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We propose a linear-optical implementation of a hyperentanglement-assisted quantum error- 
correcting code. The code is hyperentanglement-assisted because the shared entanglement resource 
is a photonic state hyperentangled in polarization and orbital angular momentum. It is possible to 
encode, decode, and diagnose channel errors using linear-optical techniques. The code corrects for 
polarization "flip" errors and is thus suitable only for a proof-of-principle experiment. The encoding 
and decoding circuits use a Knill-Laflamme-Milburn-like scheme for transforming polarization and 
orbital angular momentum photonic qubits. A numerical optimization algorithm finds a unit-fidelity 
encoding circuit that requires only three ancilla modes and has success probability equal to 0.0097. 

PACS numbers: 03.67.Hk, 03.67.Pp, 42.50.Ex 
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I. INTRODUCTION 

Quantum error correction plays an active role in the 
future realization of a quantum communication system 
[TJ [2]. Several optical experiments have already imple- 
mented simple quantum error correction routines [3 , 4 , 5J . 

The entanglement-assisted stabilizer formalism is a re- 
cent extension of the theory of quantum error correc- 
tion that incorporates entanglement shared between a 
sender and receiver [HI [7]. A further extension of this 
theory incorporates gauge qubits [8] and others give a 
structure appropriate for a stream of quantum informa- 
tion ini Uni IH]- The Ukely candidate for implement- 
ing an entanglement-assisted code is photonics because 
the entanglement-assisted model is more appropriate for 
quantum communication rather than quantum comput- 
ing. 

In this article, we propose a linear-optical implementa- 
tion of a hyperentanglement-assisted quantum code. Our 
code is hyperentanglement-assisted because it exploits hy- 
perentanglement of two photons |12j . Two photons are 
hyperentangled if they have entanglement in multiple de- 
grees of freedom such as polarization and orbital angular 
momentum (0AM) [13 E]. The benefit of hyperentan- 
glement is that a linear-optical setup suffices to perform a 
complete Bell-state analysis [THIIISIIIT]. Our proposal for 
the hyperentanglement-assisted code relies on the recent 
optical realization [TS^ of the superdense coding proto- 
col [19J and the close connection between entanglement- 
assisted quantum error correction and superdense coding 
[7]. We also employ a recent numerical optimization al- 
gorithm [20] to find an encoding circuit and a decoding 
circuit that has unit fidelity, success probability equal to 
0.0097, and requires only three ancilla modes. The cir- 



cuits act on both the polarization and 0AM degrees of 
freedom of the photonic qubits. 

We structure this article as follows. The first sec- 
tion reviews hyperentangled states, the single-photon 
polarization-OAM states, and mentions that it is possible 
to distinguish the single-photon polarization-OAM states 
with linear optics. We then discuss the superdense coding 
protocol for hyperentangled states and highlight the con- 
nection between superdense coding and entanglement- 
assisted quantum error correction. We give a brief de- 
scription of our code, its error analysis, and corrective 
operations. The final part of this article discusses the 
numerical optimization technique for finding our code's 
encoding circuit and decoding circuit. 



II. HYPERENTANGLED STATES 



The standard hyperentangled state is a state of two 
photons simultaneously entangled in polarization and 
0AM: 



- {\HH) + \VV)) ® (|00) + |00)) . 



The symbols H and V represent horizontal and vertical 
polarization respectively and the symbols O and O repre- 
sent paraxial Laguerre-Gauss spatial modes with -\-h and 
~h respective units of 0AM [^T. Changing the polariza- 
tion degree of freedom of Alice's photon in the above state 
according to the four standard Pauli operators, while 
leaving the 0AM degree of freedom unchanged, gives the 
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following four hyperent angled states: 

|$±) = i {\HH) ± \VV)) ® (|00) + |00)) , 

|*±) = ]- {\HV) ± \VH)) ® (|00) + |00)) . (1) 

We can rewrite the above four states in terms of the 
single-photon polarization-OAM states (jr^ , -0^ [30] : 

|*±) = - (±0+ (» 0± T 0A <» '/'b ± V-l ® V'b =F V^:^ «) V's) 

where A and i? label the first and second respective pho- 
tons and the single-photon polarization-OAM states 0* 
and "ip^ are as follows: 

.±_\HO)±\V O) ,±_\HO)±\V O) 



V2 



r 



a/2 



The above "quad-rail" basis states \H O), \H O), \V O), 
\V O) are four- mode single-photon states defined in 
terms of the Fock basis: 

\VO) = \Q)ho\Q)hoWvo\^)vo^ 
|yO)^|0)^^|0)^„|0)^,^|l)^,^. (2) 

Hyperentangled states are useful because a linear- optical 
analyzer distinguishes the single-photon polarization- 
OAM states <jr^ , ip^ and thus distinguishes the hyper- 
entangled states 1$='=) , 1*='=) as well il5 til6.il7] . 



III. SUPERDENSE CODING AND 

ENTANGLEMENT-ASSISTED QUANTUM 
ERROR CORRECTION 

We briefiy describe the superdense coding protocol for 
polarization encoding and hyperentangled states [18j . A 
sender Alice and a receiver Bob share a hyperentangled 
state I <&''") . Alice encodes one of four classical mes- 
sages (two classical bits) by applying one of four trans- 
formations to her half of 1$+) : (1) the identity, (2) 
\V) ^-\V), (3) \H) ^ \V), or (4) both \V) ^ - \V) 
and I if) ^^ \V). Let Z denote the second operation 
and let X denote the third operation. The result is to 
transform the original state \^'^) to one of the four 
s±\ i>Tf±\ g}jg then sends her half of the 
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states l^^" 

encoded |<I>+)'^^ to Bob. Bob performs a single-photon 
polarization-OAM state analysis in the basis (j)^,ip^ on 
each of the systems A and B to distinguish the message 
Alice transmitted. 

In the above analysis, it is important to stress that the 
dense coding transformations affect only the polarization 
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FIG. 1: (Color online) The operation of the 

hyperentanglement-assisted quantum code. Red qubits 
(those labeled as A and A\) belong to the sender Alice 
and blue qubits (the one labeled as B\) belong to the 
receiver Bob (though all qubits belong to Bob after the 
noisy channel). Alice and Bob share a hyperentangled state 
1$"'') ^ ^ prior to quantum communication. Alice uses 
the hyperentangled state to aid in encoding an information 
photon in the state \>p) . Her encoding circuit consists 
of one controUed-phase gate. She sends her photons over 
a noisy polarization-error quantum channel. Bob receives 
the photons, performs a decoding circuit, and performs two 
single-photon polarization-OAM analyses in the basis <^ , tp 
on the systems A\ and B\ to determine the error syndrome. 
Bob finally performs a recovery operation to obtain the 
information photon I?/;) that Alice first sent. 



A H 

degree of freedom of the hyperentangled state |$+) 
The classical information resides in a four-dimensional 
subspace of the 16-dimensional Hilbert space. The extra 
dimensions help in single-photon polarization-OAM state 
analysis in order to distinguish the classical messages. 

Ref. f7l discusses the close relationship between su- 
perdense coding and entanglement-assisted quantum er- 
ror correction. In superdense coding, one exploits the 
classical bits encoded in a Bell state so that Alice can 
transmit classical information to Bob. In entanglement- 
assisted error correction, one exploits the encoded clas- 
sical bits for use as error syndromes. Another way of 
thinking about this latter scenario is that Eve, the en- 
vironment, is superdense coding messages (errors) into 
the Bell states. Bob can determine the errors that Eve 
introduces by measuring in the Bell basis. 



IV. OPERATION OF THE CODE 

The operation of our code begins with an initial, un- 
encoded state consisting of one information photon and 
one hyperentangled state: 



|V')^|$^ 



AiBi 



(3) 



The information photon is as follows: 



a\H) 



f3\V) . The sender Alice possesses photons A and 
Ai and the receiver Bob possesses photon Bi. The 
entanglement-assisted communication paradigm assumes 
that Alice and Bob share the hyperentangled state prior 
to quantum communication. Figure [T] highlights the op- 
eration of our hyperentanglement-assisted quantum code. 



The sender Alice applies an encoding circuit consisting 
of one controUed-sign gate (we discuss this gate later in 
more detail) so that the state shared between Alice and 
Bob is the following unnormalized encoded state: 

She sends her photons A and Ai over a noisy 
polarization-error channel (discussed below). Bob re- 
ceives the photons and we relabel them as B and B'l 
respectively. For now, suppose that the channel does 
not introduce an error. Bob finally applies the decoding 
circuit (same as the encoding circuit) and the resulting 
decoded state is as follows: 



|$+) 



B[Bi 



(5) 



where the information photon appears in Bob's system 
B. Bob performs two single-photon polarization-OAM 
analyses in the basis (j)^ , tp^ on the systems B'l and Bi 
to diagnose the channel error. The polarization-OAM 
analysis distinguishes the four states {|<&*) , |^*)}- Bob 
measures the result !$"*") when the channel does not in- 
troduce noise. The state |$^) is a syndrome that deter- 
mines the channel error. Bob does not need to perform 
a recovery operation in this case. 



V. ERROR ANALYSIS 

In general, the channel introduces errors on the pho- 
tons that Alice transmits. We assume in this arti- 
cle that the channel is a noisy polarization-error chan- 
nel (analogous to the classical bit-flip channel.) A noisy 
polarization-error channel independently applies a polar- 
ization error X that flips the horizontal and vertical po- 
larization bases. We assume that this channel affects 
the polarization degree of freedom only and does not af- 
fect the 0AM degree of freedom. Although this channel 
may not be entirely realistic, it provides a platform for 
a proof-of-principle demonstration of the operation of an 
entanglement-assisted quantum code [5T] . 

The code protects against a single polarization error 
on either of the two photons A or Ai that Alice sends. 
It also protects against a double polarization error on 
both photons. Suppose that a polarization error occurs 
on photon A. After Bob appHes the decoding circuit, 

the state becomes X \tp) |<i>^) ^ \ Bob measures the 
photons B[ and Bi, determines they are in the state 
|<i>~), and flips the polarization of photon B to recover 
the initial information photon j?/'). Table H summarizes 
the other cases. 
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TABLE I: The table details the results of Bob's Bell state 
analysis. The states in the third column are syndromes that 
determine the channel error ("Error") and the recovery op- 
eration ("Recovery") that Bob should perform to recover the 
initial information photon. 



VI. OPTICAL ENCODING AND DECODING 
CIRCUIT 

The seminal paper of Knill, Laflamme, and Milburn 
(KLM) showed how to perform two-qubit interactions 
with linear-optical devices [22 [23]. Their method is 
a measurement-assisted scheme: it first mixes a set of 
"computational" modes and ancilla modes in a linear- 
optical device and then counts the photons in the an- 
cilla modes. The optical transformation acts on the 
computational modes. The ancilla modes help perform 
this transformation and we measure them at the end of 
the measurement-assisted scheme. The gates exploit the 
Hong-Ou-Mandel quantum interference of indistinguish- 
able photons j24]. These measurement-assisted trans- 
formations are heralded, non-deterministic, and non- 
destructive. A destructive gate involves a measurement 
on the computational modes — the informational state 
collapses to one of the states in the measurement ba- 
sis even when the gate succeeds. A non-destructive gate 
requires a measurement only on the ancilla modes — the 
result is that all the information encoded in the state 
remains intact when the gate succeeds [531 US] ■ 

The measurement-assisted scheme is useful for dual- 
rail encoded qubits and even polarization-encoded qubits 
|26j . but until now, no one has considered its appli- 
cation to "quad-rail" encoded quantum information in 
polarization-OAM states. 

The implementation of a polarization-OAM 
measurement-assisted scheme requires unconven- 
tional linear-optical elements analogous to beamsplitters 
and other tools of linear optics acting on 0AM states of 
photons. Holographic elements suffice for this purpose 
because they act on 0AM components [18j — similarly 
to the action of polarization beamsplitters on photon 
polarization. The extension of the measurement-assisted 
scheme to 0AM states is a generalization of the idea in 
Ref. [26 . There, the authors extended the measurement- 
assisted scheme to polarization-encoded qubits. Here, 
we use a similar idea to extend the measurement-assisted 
scheme to 0AM states. 

The encoding transformation corresponding to our 
code generates the encoded state in (|4|. It is a controUed- 
sign gate that acts on the four-dimensional Hilbert space 



Ha ® Hai of the information photon A and the polar- 
ization subspace of Alice's part Ai of the hyperentangled 
state in ([S]). The gate acts on the polarization degrees of 
freedom and leaves the 0AM degrees of freedom invari- 
ant. It is a linear-optical transformation on a set of six 
modes — two modes for the A system and four modes for 
the Ai system. The gate leaves the following basis states 



\Hf\H0f\ \Hf\H0f\ \Hf\V0f' 
\Hf\V0f\ \vf\H0f\ \vf\H0f\ 



(6) 



invariant and adds a phase to the remaining basis states: 



\^i 



\vr\vor' -^-\vr\v(3) 
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Ai 



(7) 



We make a statement about the mathematical struc- 
ture of the Hilbert space of polarization-OAM states. It 
is possible to decompose any Hilbert space with a tensor 
product structure. E.g., 

span{|ffO),|yO),|i/0>,|VO)} 

= spa.n{\H) , \V)} ® span {|0) , |0)} . 

While in the paraxial approximation, local optical 
transformations on the subspaces spa,n{\H) ,\V)} and 
spanJIO) , |0)} respect the tensor-product decomposition 
of the full four-dimensional space. However, the tensor- 
product notation in (ll]) from Ref. 18J may be some- 
what misleading in the context of a measurement-assisted 
transformation (which the authors of Ref. [18 do not 
discuss). A qubit-coupling measurement-assisted trans- 
formation, based on the mixing of creation operators in 
separate modes, does not respect such a tensor-product 
decomposition in general. Instead such an operation acts 
naturally on a space constructed as a direct sum, e.g., 
span{|iJ O) , 1^ O)} © spa.n{\H O) , \V O)}. The above 
restriction places a constraint on the form of the linear- 
optical encoding circuit and decoding circuit. 

Knill devised an optimal solution for the controUed- 
sign gate 157]. We could use Knill's two dual-rail qubit 
gate for an implementation of the transformation in ([7|. 
It requires the combination of two separate transforms: 
the first acts on the O 0AM states and the second 
acts on the O 0AM states. Each transformation re- 
quires two ancilla modes and has a success probability 
of 2/27 [27]. This "Knill combination" scheme thus re- 
quires four ancilla modes with a success probability of 
(2/27)^ « 0.0055. 

We have employed the numerical optimization tech- 
nique in Ref. [2U] for our encoding circuit rather than 
the above Knill combination scheme. Our optical scheme 
for the encoding transformation in (|6| and ([7| requires 
only three ancilla photons and has a success probabil- 
ity of 0.0097 (nearly a two-fold increase over the Knill 
combination scheme). 

We briefly describe the numerical optimization tech- 
nique for finding a general linear-optical circuit |20) . An 
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FIG. 2: The figure displays the results of the numerical 
gate optimization algorithm for our case of six computational 
modes and three ancilla modes. The algorithm starts an op- 
timization at a randomly selected point in the space of 6 x 6 
matrices and maximizes the gate fidelity. The second stage 
of the optimization algorithm is a constrained optimization 
of success probability within the unit fidelity subspace. The 
figure displays ordered success probability for the best 500 
optimization cycles of the total amount of 2000 optimizations 
for fidelity > 1 — 10^^. The algorithm finds an encoding cir- 
cuit with a success probability almost twice that of the Knill 
combination scheme. 



N X N matrix, where N is the number of optical modes, 
completely characterizes an optical transformation. The 
numerical implementation of the optimization algorithm 
performs a gradient search in the space of these matrices. 
The algorithm first finds a matrix that guarantees a unit 
transformation fidelity and then performs a constrained 
optimization of success probability in the unit fidelity 
subspace. Each optimization cycle ends at some local 
maximum of success probabilitiy. The resulting Kraus 
transformation, or contraction map, acting on the com- 
putational modes matches the desired target transforma- 
tion in the case of a successful measurement outcome on 
the ancilla modes. 

We now describe the results of the procedure for our 
case where the encoding circuit acts on the modes for Al- 
ice's systems A and Ai and three ancilla modes. We per- 
formed the optimization in the 81-dimensional complex 
space of a GL(9) matrix (six modes for Alice's systems A 
and Ai and three ancilla modes). Such a matrix admits a 
decomposition into a sequence of operations where each 
operation corresponds to a standard optical element [2S] . 
We were able to simplify the transformation even more 
by using a technique from Ref. |27j. The optimal trans- 
formation actually acts on three modes only: the vertical 
polarization of Alice's photon A, the vertical polarization 
and O 0AM of Alice's photon Ai, and the vertical polar- 



ization and O 0AM of Alice's photon Ai. The reduction 
of the optimal solution to a three-mode operation should 
make it easier to implement experimentally. The follow- 
ing equations illustrate the reduced transformation: 






where ni and 712 are photon numbers equal to either zero 
or one and the function / is equal to one if ni = 71,2 = 1 
and is equal to zero otherwise. One can verify that the 
implementation of the above reduced transformation is 
equivalent to the full transformation in nh. Figure ^ 
illustrates a sample distribution of success probabilities 
in increasing order. 

The best found solution provides a success probabil- 
ity of 0.00974276 at a fidelity of 1 - (6 x 10"^). This 
solution, most likely, is the optimal global solution, or 
at least close to the global optimum. One cannot verify 
the global character of a solution by numerical tools only. 
However, we have found that further optimization with 
the current scheme does not improve the result. The 
success probability shows a slight increase for non-unit 
fidelity: at the level of 0.99 fidelity, we have found a so- 
lution with success probability of 0.011. Further relaxing 
the unit fidelity requirement does not lead to any signif- 
icant increase of the success probability. 

We now address the issue of implementing our 
hyperentanglement-assisted quantum code. As described 
above, our encoding circuit solution corresponds to some 
optical transformation matrix. The first question arising 
in connection with the possible experimental implemen- 
tation of such an optical transformation is whether the 
obtained matrix is unitary. Currently, we use an opti- 
mization algorithm that searches for an arbitrary ma- 
trix, one not necessarily restricted to the subspace of 
unitary matrices. Therefore, one may need to apply a 
unitary dilation procedure j27l I29j to embed the matrix 
into a larger unitary matrix. However, a solution cor- 
responding to a maximal success probability is a matrix 
requiring minimal dilation, due to the singular behavior 
of the gradient of the success probability function on the 
manifold of unitary matrices. Indeed the two best solu- 
tions (see insert in Figure 2) correspond to matrices with 
the following singular values {l.,l.,l.,l.,l.,0.5}. One 
needs to introduce only one additional vacuum mode to 
embed the matrix into an SU(7) matrix in order to im- 
plement the transformation in the form of beamsplitters 
and phase shifters. Ref. pS; suggests the scheme of such a 
decomposition. Mathematically, it corresponds to a fac- 
torization of an arbitrary SU( A^) matrix into a product of 
SU(2) matrices. The implementation in Ref. [28' requires 
N{N — l)/2 beam splitters. However, it is well known 
that in some important cases one can find a simpler de- 
composition [27j . Knill's scheme requires only four beam- 
splitters for the realization of an SU(4) matrix, whereas 
the original KLM scheme requires eight beam splitters 
for an SU(6) matrix. The matrix we have obtained is 



rather complicated because we were not able to find a 
"nice" decomposition (See the Appendix for the matrix) . 



VII. CONCLUSION 

We have presented an optical implementation of an 
entanglement-assisted quantum code that should be re- 
alizable with current technology. The code encodes one 
information photon with the help of a hyperentangled 
state. To our knowledge, this proposal is the first sugges- 
tion for an implementation of an entanglement-assisted 
quantum code. 
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APPENDIX A 

Below we list the matrix that implements the trans- 
formation for the encoding circuit. We do not provide 
it as one large matrix because the individual entries are 
too large, but instead provide it a few columns at a time. 
The first three columns are as follows: 



-0.253936 + iO. 215424 -0.0269989 -I- i0.211134 

1 

0.0473299 + i0.183042 -0.136174 + iO.454254 



0.196523 -i0.216478 -0.233841 + iO. 184769 



-0.33549 - iO. 135251 0.314695 -I- iO. 192451 

0.318659 + i0.380869 0.3053 + ^0.314815 

0.277613 - i0.411775 -0.0145173 + ^0.484746 



The second three columns are as follows: 

-0.249991 - iO. 213976 



0.262553 + iO. 141112 

a 

-0.159651 + ^0.187183 

0a 

-0.515822 + ^0.243508 

0.220057 - iO. 296955 

0.039117 + iO. 303606 



where a = —0.611421 — lO. 791452. The next two columns 
are as follows: 



0.410744 + iO. 0245062 0.367852 - i0.184455 



-0.4806 - iO.326223 0.298488 - i0.221226 



-0.264695 - iO. 0518749 0.153492 + iO.498569 



-0.23816 - iO. 143929 -0.278242 - i0.00531807 

0.132686 - i0.193403 -0.0860369 + i0.41503 

0.382029 + i0.183596 -0.229701 + iO. 0475146 

The last column is as follows: 



-0.0349526 + iO.229345 


-0.0337073 + i0.290301 


0.403926 - iO. 202786 


-0.337872 - i0.21804 
-0.334824 - i0.268026 
-0.164744 + ^0.395987 
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